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Abstract —We consider the quantized consensus problem on 
undirected time-varying connected graphs with n nodes, and 
devise a protocol with fast convergence time to the set of 
consensus points. Specifically, we show that when the edges of 
each network in a sequence of connected time-varying networks 
are activated based on Poisson processes with Metropolis rates, 
the expected convergence time to the set of consensus points 
is at most 0{n^ log^ n), where each node performs a constant 
number of updates per unit time. 

I. Introduction 

There has been much recent interest in the design of 
control protocols for distributed systems, motivated in part 
by the need to develop protocols for networks of autonomous 
agents characterized by the lack of centralized information 
and time-varying connectivity. A canonical problem within 
the field is the so-called average consensus problem, wherein 
a group of agents must agree on the average of their initial 
values while interacting with neighbors in a (possibly time- 
varying) network. Protocols for consensus problems must be 
distributed, relying only on local information at each node, 
and robust to unexpected changes in topology. 

It is well understood by now that protocols for the 
consensus play an important role in a number of more 
sophisticated multi-agent tasks. We mention distributed op¬ 
timization, coverage control, formation control, cooperative 
statistical inference, power control, load balancing, epidemic 
routing, as examples of control and coordination problems 
with proposed solutions relying crucially on consensus. 

Our work is motivated by the observation that often 
working with real-valued variables in multi-agent control 
is neither necessary nor efficient. Indeed, limited memory 
and storage at each agent often forces the variables kept by 
each agent to lie in a discrete set. We therefore consider the 
quantized consensus problem, a variation on the consensus 
problem where the values of each agent are constrained to 
be integers lying within a certain range. Previous literature 
on this problem includes [1], [2], [3], [4], [5], [6]. As 
some of the real-world applications of the problem under 
consideration, one can consider information fusion in sensor 
networks [7], e. g., when every sensor in a sensor network 
has a different measurement of the temperature and each 
sensor’s goal is to compute the average temperature, and 
load balancing in processor networks which has various 
applications in computer science. 

The original paper [1] contained upper bounds for a 
natural quantized consensus protocol on a variety of common 
graphs. A few years later, the work of [3] proposed a 


quantized consensus protocol with upper bound of 0(n^) 
on the expected convergence time for any fixed graph. For 
dynamic graphs, [3] obtained a convergence time scaling of 
0(n®). The paper [6] obtained an upper bound of O(n^), 
but only on complete graphs. Recently, it was shown in 
[8] that a certain “unbiased” quantized consensus protocol 
has maximum expected convergence time of 0{mnD logn) 
on static networks where m and D denote the number 
of edges and the diameter of the network, and has maxi¬ 
mum expected convergence time 0(nTOmax-Dmax log^ n) on 
connected time-varying networks, where Wmax and Z^max 
denote, respectively, the maximum number of edges and 
the maximum diameter in the sequence of time-varying 
networks. Unbiasedness means that each step of the protocol 
was based on choosing a random edge uniformly among all 
possible edges in the network. 

A faster upper bound on convergence time and only for 
static networks was given in the recent paper [9], where 
a protocol was provided whose expected convergence time 
in general static networks is O(n^logn) Q As of the time 
of writing this paper, the upper bounds of 0(n^ logn) and 
0(nmmaxT^max log^ n), respectively, are the fastest proto¬ 
cols known to us for randomized quantized consensus over 
static and dynamic networks [8], [9], respectively. It is worth 
noting that the convergence speed of a protocol is measured 
by the maximum over all initial inputs of the expected time 
that the given protocol will run until reaching consensus, 
where each node performs a constant number of updates per 
unit time. 

In this paper, we analyze a protocol wherein nodes coop¬ 
erate to perform updates on edges connecting them at so- 
called Metropolis rates. We are able to show that when each 
node performs a constant number of updates per unit time, 
convergence time on connected time-varying networks is at 
most 0{n? log^ n), giving us the fastest quantized consensus 
convergence guarantee over time-varying networks. 

The paper is organized as follows. In Section [I^ we 
formulate the quantized consensus problem and discuss some 
of its properties. In Section III we show that quantized 
consensus with Metropolis rates has essentially quadratic 
convergence time on fixed networks; this result is not better 
than the previous convergence time of [9], but we include 


it here because our later results rely on it. In Section IV 


we prove our main result, namely that Metropolis quantized 


*The bound given in [9] was 0(n^ logn), but this was a count of the 
total number of updates; in terms of time, this leads to a quadratic expected 
convergence time. 




consensus converges in essentially quadratic time on time- 
varying networks. We conclude the paper by identifying 
some future directions of research in Section FV] 

Notation; We let [n] = {1,2,... ,n}. For an undirected 
graph Q = (V, £) we let N{x) be the set of neighbors of x, 
and furthermore, let d{x) denote the degree of x i.e., d{x) = 
|iV(a;) |. For a given vector v, we denote its *th entry by Vi and 
its transpose by v'. We say that a matrix A with nonnegative 
entries is stochastic if each row of A sums to 1. We say A 
is doubly-stochastic if both A and A' are stochastic. 

II. Problem Formulation 


of an activated edge have integer values which differ by at 
least 2. 

It has been shown earlier in [1] that such dynamics 
based on update rule 0 will converge with probability 1 
to consensus set defined by 

C = {x|T,e{Lx(0)J, LS(0)J+1}, = x(0),z =€ [n]}. 

In the remainder of this paper our goal is to provide an upper 
bound on the expected convergence time of the quantized 
Metropolis dynamics to the consensus set. 


In this section, we assume that we are given a fixed, 
undirected, connected graph Q — {V,£) with \V\ = n and 
without self-loops. The Metropolis Markov chain on this 
graph is then defined as follows. 

Let be a square matrix whose ij’th entry is defined as 

f max{d{i) ,d{j)} ’ if j, j e N{i) 

~ I ~ if * = J (f) 

[O otherwise, 

Note that Ai is symmetric, nonnegative, and doubly stochas¬ 
tic. We refer to the Markov chain which transitions according 
to Ai as the Metropolis chain. Moreover, given nodes x,y G 
V, the hitting time H'^{x,y) is the expected time until the 
Metropolis chain with initial position at node x reaches node 
y (quantities associated with the Metropolis chain will gener¬ 
ally be denoted with an “m” superscript). For future use and 
convenience, we adopt the notation A^- = 1/ max((i(z), d{j)) 
whenever {i,j} G £ and = 0 otherwise. 

A. Quantized Metropolis dynamics 

We next introduce a continuous time quantized process 
based on Metropolis weights, whose behavior will turn out 
to be related to the Metropolis chain. For each link {i,j} G £ 
of the graph Q, we consider a Poisson process of rate A^ . 
Each node i begins with an integer value Xi(0) which lies 
in the range \l,L]. Each time the process corresponding to 
an edge {i,j{ G £ registers an arrival, the two nodes i,j 
perform the quantized consensus update from [1]: 

( Xiit) - 1, if x^(t) > Xj{t) 

Xi(t) + 1, if Xj(t) < Xj{t) (2) 

Xi{t), if x,{t) = Xj{t), 

and likewise for node j. In other words, each Xi{t) is an 
integer-valued jump process whose jumps occur whenever 
an arrival occurs at any edge incident on i. 

Note that the above update rule allows for the possibility 
that Xi = Xi{t). In this case, we will say that the update 
at time t was trivial. Eurthermore, observe that if \xi{t) — 
Xj{t)\ — 1, then the update of Eq. (|^ will cause nodes i 
and j to swap values, i.e., Xi{t~^) = Xj{t),Xj(t'^) = Xi{t). 
In this case, we will also say that the update was trivial. If 
neither of these two cases has occurred during an update, we 
will say that the update was non-trivial. Simply speaking, a 
non-trivial update refers to the case where the incident nodes 


III. Convergence Time over Static Networks 

We now begin the analysis of the convergence time of 
the quantized Metropolis algorithm over a fixed graph Q. 
The bounds derived in this section are of the same order of 
magnitude as the previous work [9], and the proof follows 
similar lines as well; however we are including it here as 
our analysis over time-varying networks in the next section 
relies on these results. 

A key step in our analysis is to bound the expected 
time until the first nontrivial update takes place. Eor this 
purpose, we let T be the maximum expected time such that 
a nontrivial update takes place over all possible configura¬ 
tion of integers. More precisely, for integers ji,... ,jn let 
T{ji,... ,jn) be the expected time until a nontrivial update 
takes place when node i begins with integer value ji. Then 


T = max 




Definition 3.1: Consider two random walkers moving 
based on whether the activated edge in the quantized 
Metropolis dynamics is incident to them. That is, if one of 
the walkers is at node x, and if the next edge to register an 
arrival is incident to x, i.e., if it is {x, y} for some y G N{x), 
then the random walker moves from x to y. We refer to such 
a process as the original process. We denote by AI°{x, y) the 
expected “meeting time” of this process, defined to be the 
expected time until an edge incident to both walkers registers 
an arrival provided the two walkers started at nodes x, y with 
X 7 ^ 2 / (in general, we will denote quantities associated with 
this process with an “o” superscript). By convention, we set 
M°{x, x) = 0 for all x gV. 

Proposition 3.2: T = m.a'x.x^y^v M°{x,y). 

Proof: The proof is immediate and is not included due 
to space limitation. ■ 

Based on the above proposition, our next step is to bound 
maxx^y^v M°{x, y). We will actually find it easier to instead 
bound a meeting time associated with a slightly different 
process, which we call the virtual process, defined next. 

Definition 3.3: The virtual process is identical to the orig¬ 
inal process until the two walkers x, y become each other’s 
neighbors in Q. At that time, the edge connecting them 
registers arrivals according to a Poisson process of rate 2Axy. 
We denote the expected meeting time function of the virtual 
process by M''{x,y) (in general, we will denote quantities 




associated with this process with a “v” superscript). By 
convention we set x) = 0 for all x gV. 

Definition 3.4: A function ft, : ^ K is called harmonic 

at a vertex x G for a Markov chain with transition 
probability matrix P if ft(x) = '^y^^P{x,y)h{y). 

Remark 3.5: Given a nonempty subset B C D. and a 
Markov chain with an irreducible transition matrix P, every 
harmonic function ft(-) : —>■ K over \ B which satisfies 
ft(x) > 0, Vx G must be nonnegative over the entire (see 
Proposition 9.1 of [10]). 

Definition 3.6: A vertex 6 is called a hidden vertex of the 
Metropolis chain if H^{0,x) < ift“(x,6*),Vx G V. 

Remark 3.7: It is known that the hitting times of a re¬ 
versible Markov chain satisfy transitivity property, i.e., 


H^{x, y) z)+H^{z, x) = H^{x, z)+H^{z, y) +H^\y, x), 


see [11] for a proof. As a consequence every reversible 
Markov chain has at least one hidden vertex [12]. 

In the following lemma we bound the expected meeting 
time of the virtual process by showing that the expected 
meeting time function of the virtual process M''{x,y) satis¬ 
fies essentially the same recursion as the function $“(•, •) : 
y X IG —K defined by 

$“(x, y) = ift“(x, y) + 9) - y), (3) 


where ft is a fixed hidden vertex of the Metropolis chain. 
Lemma 3.8: For all x,y GV, we have M''{x,y) < 6n^. 
Proof: Fix nodes x and y such that x f y, and define 


^ ^ ^xxi: ^ ^ ^yVi i ^xy — ^x 


- A„ 


Xi£N{x) 


yi&N(y) 


An immediate consequence of Remark 13.71 is the symmetry 
of i.e., u), Vu, v G V. Therefore 

we can argue that 


+ 


^x\ 

-A..T? 


.4\.x 

xi^N(x) 




yiSN(y) 


A. 


+ ^[H-iy,0)-H-ie,y)] 


A 


+ ^[77”(x,ft)-7f”(ft,x)], 


A 


where in the second equality we expanded H'^{x,y) and 
ift“(y,x). Simplifying the last relation, we obtain 


<I-”(x,y) = ^+ ^ 

XieN(x) ‘^^y ViSNiy) ^'-^y 

+ ^ 0) y)) + ^ e)-H-^{e, x)), 

i^xy 4\xy 


(4) 

Using the definition of $“(•,•), for each Xi G N{x) 
we have, <i)™(xi,y) = H'^{xi,y) +iT“(y, ft) — ift“(ft,y). 
Multiplying this relation by and summing over Xi G 

4\.xy 


N{x), we obtain 

xi&N{x) x,&N(x) ^^^y 

+ ^[iJ“(y,ft)-iT“(ft,y)]. (5) 
By symmetry of and using the same argument we get 

^ ^ci>“(x,y,)= E 

y,GNiy) '^-y y,GNiy) 

+ ^[ift“(x,ft)-ift“(ft,x)]. (6) 


xy 


Substituting and Q in 0 yields 

XiSN^x) 


VieN(y) 


(7) 


On the other hand, we note that the meeting time of the 
virtual process can be expanded as 

M''{x,y) = ^+Y ^Af''(xi,y)+ ^ ^^M''{x,yi). 

^^^y x^^Ntxf^^y yi€N(yf^^y 

(8) 

We therefore see that i$“(x,y) and M^{x,y) follow 
the same recursion formula when x f y. Thus defining 
/(x, y) = ^$“(x, y) — M''(x, y) we have 'ix f y gV 


f{x,y)= E E 

x,eNix) ^y 


v.^Niv) ^^y 


Defining the stochastic irreducible matrix Q as 

ifs = y, rGN{x) 
Q((a;,y), (as)) = S Ei’ ifr = x,wGN{y) 

[o, Else, 

one can see that f{x,y) is harmonic for the matrix Q over 
U X U \ {(x,x) I X G U} and, since ft is a hidden vertex, 
we have that for all x G V, /(x,x) = ^(j 7™(x,ft) — 
ift“(ft,x)) > 0. Therefore, using Remark 3.5 we imme¬ 
diately get f{x,y) > 0 for all x,y G V, and therefore 
My{x,y) < i$™(x,y). Finally, it was shown in [13] that 
H'^{x,y) < Qnf for all x,y G V, which now immediately 
implies the current lemma. ■ 

Our next lemma shows that max^, j, M°(x, y) is essentially 
of the same order as max^, ^ M''(x, y). 

Lemma 3.9: uiaxx^y M°{x,y) < 12n^. 

Proof: Fix x f y and initialize both the original 
process and the virtual process with initial positions of the 
walkers being x and y. Furthermore, let us couple these 
processes to move identically until the first time when the 
walkers are each other’s neighbors. At this time assuming 
that the walkers are at nodes u and w, let us split the edge 
connecting them of rate 2A„u, in the virtual process into 
two distinct edges of rate (Figure [^. We will make 
the first of these (the bottom black edge in the right side 
of Figure [T]), as well as the edges going from u and w to 
their neighbors, register arrivals in the virtual process exactly 









when the corresponding edges in the original process register 
arrivals. Now, for these coupled processes, let us denote by 



Fig. 1. Interpretation of the virtual process using the original process and 
an extra edge. 


T^{x,y) a random variable denoting the first time when the 
random walkers in the virtual processes meet and T°{x,y) 
the first time when the random walkers in the original process 
meet, given that the processes have been initialized from x 
and y. We have that T°{x,y) = T^{x,y) if the first edge 
of rate X^w registers an arrival before the second edge of 
rate Xuw This happens with probability If this does not 
happen, T°{x,y) is T''{x^y) plus another random variable 
which is a meeting time in the original process starting from 
u and w, i.e., t°{u,w). In other words. 


T°ix,y) = 


'{x,y), 


w.p. i (black edge ticks) 


T^{x,y) + T°{u,w) w.p. i(red edge ticks). 
Taking expectation of the above relation, we obtain 

M°{x,y) < l^M'^ix.y) + i(M''(a;,y) + maxM°(u,ui)) 

Z Z u,w 

which implies that max^, j, M°(x, y) < 2maXx^y M''{x,y) 
which by Lemma 3.8 is at most 12n^. ■ 

Remark 3.10: In view of the previous lemma, we have 


P(T°(a;,y) > f) < e Li2e,i2J^ 


(9) 


This can be seen by the same argument as in Chapter 2.4.3 
of [11] whose argument applies verbatim here. 

Theorem 3.11: The expected time until the quantized 
Metropolis dynamics reaches a consensus set is 0{n^ logn). 

Proof: Consider the time it takes for the Lyapunov 
function V (t) = max^ ^^(f) — min^ Xi(t) to shrink by at least 
1 starting from any non-consensus configuration. We show 
that this time is 0(n^ logn), which proves the theorem. 

Indeed, let S'inax(f) = argmax* a;i(f), Smin(t) = 
argmini Xi(t). Now let us consider the event that V(t') = 
V(t) for t' > t. We claim that there must exist i G 
Sniax{t),j G 5'min(f) such that two random walkers in 
the original process starting from nodes i and j have not 
yet met during the time interval To see this, we 

note that during the executions of the quantized Metropolis 
dynamics no new values maxiXi{t) and mmiXi{t) will 
be created; these values only travel from one node to the 
other over the network. Now given that at time f we have 
V{t') = V{t), there must exist two nodes i', j' such that 
= maxiXi{t) and Xj'{t') = mmiXi{t). By tracking 
back the origin of the random walks which bring the values 
maxi Xi{f) and min^ Xi{t) at time t to the nodes i' and f at 
time t', one can see that these walks have never meet, and 
they must have originated from some nodes i G G 


S'min(f)- Therefore, the probability that V(t') = V(t) for 
t' > t is upper bounded by the probability that there exist 
some i G Syna^{t),j G S'min(f) such that two random walkers 
in the original process starting from nodes i and j have not 
yet met. Thus letting C being the time elapsed until V{f) 
shrinks by 1, we have that 


£((7) = J F(C>t)dt< J min 1 1, ^ P(T°(a;, y) > f)| df 

POO 

< J min |l, i 2 ^)| df = O logn) , 

where the second inequality is due to Eq. 0 and the last 
equality follows from min {1, Ae”®*} df = , 

which holds when A> I, and is from [11], Ch. 5.3.2. ■ 


IV. Convergence Time over Dynamic Networks 

In this section we analyze the expected convergence time 
of the quantized Metropolis dynamics over time-varying 
networks and prove our main result, namely that this conver¬ 
gence time is essentially quadratic in the number of nodes. 

Toward this aim, let us consider a sequence of connected 
undirected networks Q{t) = {V,£{t)),t > 0, over the 
same set of vertices V which may change at discrete time 
instances, i.e., Q{t) = ^([tJ),Vf > 0. Thus our notation 
for the set of neighbors of x will now be Nx{t), which 
depends on t. In contrast to the previous section, we will 
now assume that each node always possesses a self-loop, 
i.e., X G N^{t) for all a; G V, f > 0. We thus introduce the 
notation N^{t) to denote the neighbors of node x excluding 
itself; N^{t) = Nx{t)\{x}. The degree of a node x, denoted 
by dx{t), will be referred to as the cardinality of N^{t). 

A. Quantized Metropolis Model over Time-Varying Networks 

Given a network Q{t) at time f > 0, we associate a 
Poisson process with each edge {i,j} of Q{t) with a rate 
of Xij{t) = 1/max{di{t),dj{t)), i.e., the Metropolis weight 
corresponding to that edge at time t. When an edge {i,j} G 
£{f) registers an arrival, we let the incident nodes update 
their values based on 0- Moreover, for each node x and 
time t, the self-loop {x, x) registers arrivals according to a 
Poisson process with rate of 

Xxx{t) = l-]^ Y, (10) 

XiGV'(t) 

Note that X^x (t) > 0 for all x G V, f > 0, and the sum of 
the rates of the edges of Q{t) is always equal to \V\ = n. 
When a self-loop registers an arrival, no update is made. 

B. Preliminary Definitions and Relevant Results 

We next introduce the notions of the “original” and “vir¬ 
tual” processes over the time-varying graph sequence {G{t)}. 

Definition 4.1: Consider two random walkers moving 
based on whether the activated edge in the quantized 
Metropolis location is incident to them or not. That is, if 
a random walker is at node i at time t when the edge 
{fj} G G{t) registers an arrival, the walker moves to node 











j. As before, we will refer to this as the original process, 
but note that the graph sequence is now time-varying. The 
meeting time M°{x,y) for x y is the expected time until 
two random walkers, starting at locations x and y at time t, 
are both incident on an edge which registers an arrival. By 
convention we set M°{x, x) = 0 for all x G V,t > 0. 

Definition 4.2: We dehne the virtual process to be iden¬ 
tical to the original process except when the two walkers 
in the original process are each other’s neighbors at nodes 
u and w for some t > 0, i.e., u G N.w{t)- At this time in 
the virtual process we let the edges {r, s} G S{t) register 
arrivals at rates y,rs{t) dehned as follows. 


Continuing with ( [T3| ) and since Ti is exponential with 
parameter n, we further have that 

/ oo 

= i (^1 - Ml{x,y)dU 

+ (14) 


ftrs(f) — 



if {r, s} = {u, w} 
if {r, s} = {u, u} 
if {r, s} = {tu, w} 
if else. 


( 11 ) 


Note that fj,rs (f) > 0 for all r, s S V, f > 0. We refer to 
the meeting time of the virtual process starting from x and 
y at time t as T^{x,y) and its expectation M^{x,y). By 
convention, we set {x,x) =0 for all x GV,t > 0. 

In words, the rate of the edge {rt, u>} is doubled, while 
the rate of the self-loops is decreased to make sure that all 
rates still sum up to n. 


C. Convergence Rate over Time-Varying Networks 

In this part, we state our main results for the quantized 
Metropolis dynamics over time-varying networks. Dehne 
h.xy{t) to be 

^xy {t)= E! (t) + E! ^Vi (0 =2(2 — Xxx (t) ~ Xyy (t)) . 

( 12 ) 

Fix t and let Ti be the hrst time after t that an edge registers 
an arrival in the virtual process. We then have 


E [rj(x, y)\n] = (Ti - t) + (^1 - {x, y) 


Note that this equation holds regardless of whether x and y 
are neighbors, and this fact is the reason why we introduced 
the virtual process in the hrst place. Next the following 
lemma shows that Ti follows an exponential distribution with 
parameter n. 

Lemma 4.3: Let us consider a process obtained by restart¬ 
ing a Poisson process of rate n at every discrete time instant 

k = 0,1,2,.The hrst time arrival Ti of the new process 

thus generated will follow an exponential distribution with 
parameter n. 

Proof: The proof follows by computing the cumulative 
function of Ti and showing that P(Ti < f) = 1 — 


Now let us dehne v{f) to be a column vector of length 
n{n — 1) whose entries are the variables Mf {x,y),Vx y 
in any order. It follows that we can write the above recursion 
for M'l (x, y) in matrix form as 


v{fi) = -l-f 
n 



ne D (ti) V (fii) dti, 


(15) 


where D{t) G i)xn(n i) ^ matrix whose rows and 

columns we will index by (x, y), x 7 ^ y as 





^xy(i) 


if s = y,r G Nx(t) \ {x,y} 
if r = x,s G Ny{t) \ {x,y} 
if (r,s) = (x,y) 

Else. 

(16) 


We change variables in ( fTS] ) from to z = —e 
and obtain 


.(() = -!- 


'-1 


D t- 


In(-z) 


V t — 


In(-z) 


dz. (17) 


We justify this change of variables by appealing to Theo¬ 
rem 2631 in [14]. Indeed, the equivalence of Eq. ( [T7) i and Eq. 
(15 I is an instance of the equality g = Jj g{cj){z))(l)'(z). 
Here (j){z) = —while g{z) = D{t—ln{—z)/n)v{t— 
ln(—z)/n). Theorem 2631 in [14] allows us to assert this 
equality subject to (i) g{-) being Lebesgue measurable (ii); 
and (/)(z) being absolutely continuous on any closed bounded 
subinterval of I. Item (ii) is clearly satished here. Item (i) 
follows because D{-) is a piecewise continuous function by 
dehnition, and v{-) is a continuous function (indeed, taking 
tn ^ t and conditioning on no transitions in , t] as well as 
no meeting occurring in the same interval, we immediately 
obtain the continuity of v{t) as a function of t. 


Next, we note that it is immediate from Eq. ( [Tbl l that D{t) 
is a symmetric matrix. Furthermore, it is easy to see that 
D{t) is sub-stochastic, implying that its eigenvalues are all 
real and less than 1 in modulus. Our next step is to upper 
bound both the largest and smallest eigenvalues of D{t). 

















The first step is to extend D{t) to a stochastic matrix 
P g as follows: 


1, 

if X = y = r = s 

A. At) 

n ’ 

if x/y,s = y,r GN^{t)\{x} 

Asy (t) 

n ’ 

if x/y,r = x,s GNy{t)\{y} 

1 Axy(i) 

n ’ 

if x/y, {r,s) = {x,y), 

0 , 

if else. 


(18) 


Note that since P{t) is stochastic, it can be interpreted as 
the transition matrix of a Markov chain with absorbing states 
S = {{x,x) : X G V}. When we adopt the convention that 
states in S correspond to the first n rows of P{t), we have 
that the matrix P{t) can be represented as 


P{t) = 


In X n 0 j 

C{t) D{t) ) ’ 


(19) 


where C{t) and D{t) are, respectively, matrices of sizes 
(n^ — n) X n and (n^ — n) x {n? — n). 

Lemma 4.4: Consider a Markov chain with transition ma¬ 
trix Q = ( ^ z!) ) ’ additional property 


that there is a path starting from any node and ending in the 
first k nodes. Furthermore, let us assume that D is symmetric 
and let us denote its largest eigenvalue by Ai„ax(^)- Let 
Hi denote the expected time until the chain is absorbed in 
{1,..., fc} starting at node i and H = max^ Hi. Then 


Amax(-D) <1-4- 
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Proof: Indeed, 1/(1 — Xmax{D)) L a positive number 
which is an eigenvalue of (/ — D)~^ and is consequently 
upper bounded by ||(/ — Zl)“^||oo- But by [15], Theorem 
4.5, we have that the sum of the entries in the f’th row of 
(/-£»)-t is iZ,. ■ 

Lemma 4.5: For all f > 0, Amax(-D(f)) < 1 - 

Proof: Fix t and let H^^^\{x,y}) denote the hitting 
time to S' = {{f, i} | f G F} in P{t) starting from {x, y}. It 
is immediate that when x f y, 




E 


VieN'it) 


\vi (^) 

I^xy(t) 


H^^^\{x,y.}). 


Comparing this to Eq. and noting that H^^*^{{i,i}) = 0 
for all f g y, we obtain that H^^*^ {{x,y})/n equals the 
expected meeting time in the virtual process on the graph 


G{t) starting from x and y. Appealing to Lemma 3.8 


4.4 


then 


can conclude that expected hitting time in P{t) t o the first 
n states is at most 6n^. Appealing to Lemma 
completes the proof. 

Observe that the lower bound 


Xmin{D{t)) > 1 - - 
n 


( 20 ) 


follows immediately by Gershgorin circles due to the obser¬ 
vation that Axy{t) < 2 for all x,y G V,t > 0. We are now in 
a position to state and prove the main result of this section. 

Theorem 4.6: The expected time until the quantized 
Metropolis dynamics over time-varying connected networks 
reaches the consensus set is 0(n^ log^(n)). 

Proof: Note that we may iterate the recursion of Eq. 
0 to obtain 


t— 


t - 

ln(- 2 ) ln(- 2 ') 


1 f ln(— 2 ) A / 1 

v{t) = - 1 + / D[f ' ' 

+/>( 


•V t — 


\n{-z) ln(-20 


dz' ]dz 


2 

< -1 
n 




/ U rO 

-J-i 


Dlt- 


In(-z) ln{-z') 


■ V [ t — 


ln{-z) ln(- 2 ') 


dz'dz, 


where the last step used the sub-stochasticity of D{-). 
Iterating this k times, we obtain 


u(f)< -l-b 
n 


[- 1 . 0 ]* 


j=i 




( 21 ) 


Now let us use introduce the notation M = sup^Q ||u(f)||oo- 
Elementary arguments suffice to establish that M < 00 
recall that v{t) stacks up expected meeting times in the 
virtual process starting at time t, and the finiteness of M 
can follow from the observation that the probability of 
intersection in any positive length interval [t, t+a] is positive 
and bounded away from zero independently of t. Thus for 
every e > 0 there exists t* such that M — e < ||'(;(f*)||oo- 
Since the entries of D{-) and v{-) are always nonnegative, 
we can take the infinity-norm from both sides of (|2T| and 
use triangle inequality to obtain 


M-e< ||t;(f*)||oo < - 



* v^ ln( —Zt) 


i=l 


dzk-.-dzi 


^ln(- 


< -+M 
n 


<- +M 
n 


i{n — 1) 


2 = 1 
k 


dzk-.dzi 


1=1 


n{n-l) n 
1=1 


2=1 

j 




ln(-Zi 


i=l 


dzk...dzi 


dzk-.dzi 


< - + My/ n(n - 1) ( SUp{|Amax(-D(f))|, |Aniin(L’(f))|} 


\t>0 




( 22 ) 


where in the fourth inequality we have used the fact that 
the infinity-norm of a matrix is always upper bounded by 
its induced 2-norm times the square root of its dimension, 
i.e., ||iA(-)||oo < y/nfa — l)||Zl(-)|j2. Moreover, the last 
inequality is valid due to Lemma |4.5| and Eq. 

























































Let us choose k* so that \Jn{n — 1) (l — < \- 

Appealing to the inequality (1 — IjxY < e“^, we obtain 
that we may choose k* = 0{nP \ogn) to accomplish this. 
Plugging this into the last line of Eq. we immediately 
obtain M — e < ^ Since this holds for any e > 0, it 
implies that M < 2k*/n = 0{v? logn). 

We have thus obtained an upper bound of O(n^logn) 
on the expected meeting time in the virtual process starting 
from any pair of nonidentical nodes and any time t. The 
rest of the proof directly parallels the arguments we have 
made in the fixed graph case, and we only sketch it rather 
than repeat the relevant parts verbatim. The first step is to 
argue that 2M is an upper bound on the meeting time of the 
original process. The proof proceeds exactly as in the case of 
Lemma [3^ by coupling the two processes and conditioning 
on the transition at the last step in the meeting time of the 
virtual process. The next (and final) step is to argue that the 
time it takes V{t) = maxia;i(f) — vcmvixYt) to shrink by 
1 is upper bounded by 0(n^ log^n). This follows exactly 
as in the proof of Theorem |3.11| Indeed, we can argue that 
if V{t') = V{t) at some time t' > t, this means that a 
pair of nodes performing a random walk according to the 
original process have not met between times t and t'. Upper 
bounding the latter using the union bound, we once again 
obtain that the expected time it takes for V(t) to shrink by 
1 is the expected meeting time in the virtual process times 
a multiplicative factor of O(logn). ■ 

V. Conclusion 

We have studied the quantized consensus problem on 
undirected connected networks in both static and time- 
varying settings. In particular, we have proved an upper 
bound of 0(n^ log^ n) for the convergence time of quantized 
Metropolis over time-varying connected networks. 

A future direction of research would be to improve con¬ 
vergence times further. For example, [16] attained a linear 
convergence time for consensus on any fixed graph, and it is 
an open question to obtain a quantized consensus protocol 
which replicates this. Moreover, an interesting problem is to 
see to what extent the results here can be carried over to 
protocols with nonlinear transmission [17]. 
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